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(N ■ ABSTRACT 

O . Let X be a regular linear diffusion whose state space is an open interval C M. 

I We consider a diffusion X* which probability law is obtained as a Doob h- 

O ■ transform of the law of X, where /i is a positive harmonic function for the 

. infinitesimal generator of X on E. This is the dual of X with respect to 

|Vj I h{x)m{dx) where m{dx) is the speed measure of X. Examples include the 

^ ■ case where X* is X conditioned to stay above some fixed level. We provide a 

\ construction of X* as a deterministic inversion of X, time changed with some 

, random clock. The study involves the construction of some inversions which 

^ . generalize the Euclidean inversions. Brownian motion with drift and Bessel 

>■ 

■ processes are considered in details. 
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1. Introduction and preliminaries 

Let X := {Xt,t < () he a regular diffusion with life time ( and state space E = 
(Z, r) C M which is defined on complete probability space (J^()t>o,P). Assume that X 
is killed, i.e. sent to a cemetery point A, as soon as it hits one of the boundaries; that is 
( = inf{s, Xs = I OT Xs = r} with the usual convention inf{0} = +00. 

Our objectives in this paper are summarized as follows. Given a positive function h 

which is harmonic for the infinitesimal generator L of X, i.e. Lh = 0, we give an explicit 

construction of the dual of X with respect to h{x)m{dx) where m{dx) is the speed measure 

of X. The dual, which distribution is obtained by a Doob /i-transform change of measure 

of the distribution of the original diffusion is either the process itself, i.e. the process is 

self-dual, or the original diffusion conditioned to have opposite behaviors at the boundaries 
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when started from a specific point xq in the state space; this is explained in details in 
Proposition m below. We refer to the original paper [5] by J.L. Doob for /i-transforms and 
to [3] where this topic is surveyed. The procedure consists in first constructing the inverse 
of the diffusion with respect to a point xq E E which is a deterministic involution of the 
original diffusion. Time changing then with an appropriate clock gives a realization of the 
dual. The construction is known for the case where X is a three dimensional Bessel process 
started at a positive Xq. In this case, the inverse process with respect to Xq = 1 is {1/Xt) 
and the dual process is a Brownian motion killed when it hits 0, see [H]. Of course, the 
former process is a Brownian motion conditioned via a Doob /i-transform to stay positive. 
Similarly, the three-dimensional hyperbolic Bessel process can be realized as a Brownian 
motion with negative unit drift conditioned to stay positive. In fact, this work is inspired 
by the discovery of the inversion | In coth Xt of the three dimensional hyperbolic Bessel 
process {Xt) with respect to Xq = | ln(l + -\/2); when this is appropriately time changed, 
we obtain a Brownian motion with negative unit drift. The motivation of the search for 
deterministic inversions between stochastic processes comes from Potential Theory, where 
the Kelvin inversion plays a crucial role, see e.g. [2]. 

In order to say more, let us fix the mathematical setting. Suppose that X satisfies the 
s.d.e. 



where Xq G {Wt^t < C) is a standard Brownian motion and a, 6 : — > M are 
measurable real valued functions. Assume that a and h satisfy the Engelbert-Schmidt 
conditions 



where L]^^{E) is the space of locally integrable functions on E. Condition fll.2p implies 
that (11. ip has a unique solution in law, see Proposition 5.15 in [H]. We write {J^^,t > 0) 
for the natural filtration generated by X and denote by VJ^{E) the set of diffusions 
satisfying the aforementioned conditions. For background on diffusion processes, we refer 



to ( [3], [9], mm)- 

Let X G VF{E). For y E E, let Hy = mi{t > 0; Xt = y} be the first hitting time of y 
by X. Recall that the scale function of X is any continuous strictly increasing function 






a 7^ and b/a^ G LUE), 
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on E satisfying 



:i.3) 



F,{H^ < Hp) = {s{x) - s{P))/{s{a) - 



for all/<a<x</3<r. This is a reference function which is strictly increasing and 
given modulo an affine transformation by = J^^exp (—2 h{r) / a"^ {r)dr)dz for some 
c & E. For convenience, we distinguish, as in Proposition 5.22 in [9], the following four 
different subclasses of diffusions which exhibit different forms of inversions. We say that 
X e VJ^{E) is of 

• Type 1 if — oo < s{l) and s(r) < +oo; 

• Type 2 if — oo < and s(r) = +oo; 

• Type 3 if s(/) = — oo and s(r) < +oo; 

• Type 4 if = — oo and s(r) = +oo. 

Type 4 corresponds to recurrent diffusions while Types 1-3 correspond to transient ones. 
Recall that the infinitesimal generator of X is given by Lf = ((t^/2)/" + bf where / is 
in the domain T>{X) which is appropriately defined for example in [8]. For xq & E let h 
be the unique positive harmonic function for L satisfying h{xo) = 1 and either 



h{l) 



or 



l/h{r) if X is of type 1 with 2s{xo) ^ s{l) + s(r); 

1 if X is of type 1 and 2s (xq) = s{l) + s(r) or type 4; 

if X is of type 2; 



h{r) = if X is of type 3. 



If X is of type 4 or of type 1 with 2s (xq) = s(/) + s(r) then /i = 1 otherwise h is specifically 
given by (12. 7p which is displayed in Section 3 below. 

Let X* be the dual of X, with respect to h{x)m{dx) , in the following classical sense. 
For all t > and all Borel functions / and g, we have 



f{x)Ptg{x)h{x)m{dx) = / g{x)P^ f{x)h[x)m{dx) 



where Pj and are the semigroup operators of X and X*, respectively. The probability 
law of X* is related to that of X by a Doob /i-transform, see for example |6] . To be more 
precise, assuming that Xq = x E E then the distribution P* of X* is obtained from the 
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distribution P^,. of X by the change of measure 



(1.4) 



dF: 



>* I 



dF. 



t<C- 



h{x) 



We shall denote by E* the expectation under the probability measure P*; X* has the 
infinitesimal generator L*f = L{hf)/h for / G V{X*) = {g : E ^ E, hg e V{X)}. The 
two processes {h{Xt),t < () and {l/h{X^),t < (*) are continuous local martingales. We 
shall show that the former process can be realized as the latter when time changed with 
an appropriate random clock. Thus, the expression of either of the processes X* and X, 
which are both i?-valued, in terms of the other, involves the function I{x) = h^^{l/h{x)) : 
E ^ E which is clearly an involution of E in the sense that / o /(x) = x, for all x & E, 
and I{E) = E. For general properties of real valued involutions, we refer to [20] and [21] . 
Now, indeed, a variant of our main result states that when h is not constant the processes 
{h{X^),t < (*) and {l/h{Xrt),t < A^), if both are started at xq, are equivalent, where, 
for t > 0, 



and Tt is the inverse of At. Two interesting features of the involved clocks are described 
as follows. First, Tt = where AJ' is defined as above with X replaced by X*. Second, 
we have (* (resp. () and A^ (resp. A^,) have the same distribution; these identities in 
distribution for killed diffusions resemble the Ciesielski- Taylor and Biane identities, see 
[1] and [1]. We call the process {I{Xt),t < Q the inverse of X with respect to xq. 

In the transient case, the general construction of X* from X discovered by M. Naga- 
sawa, in [12], applies to linear diffusions; see also [H] and [11]. We mention that this 
powerful method is used by M.J. Sharpe in [19] and by D. Williams in [22] for the study 
of path transformations of some diffusions. While the latter path transform involves time 
reversal from cooptional times, such as last passage time, the construction we present 
here involves only deterministic inversions and time changes with random clocks of the 
form (11. 5p . While we only consider one dimensional diffusions in this paper, inversions 
of stable processes and Brownian motion in higher dimensions are studied in [2] and [23] , 
respectively. 

This paper is organized as follows. We start section 2 with the construction of a family 
of inversions associated to a diffusion X G 'DJ^{E). The construction involves a reference 
scale function s and not the speed measure of X. We note that the inversion of E in 



(1.5) 
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the direction of s is uniquely characterized by the fixed point Xq. Also, among the set of 
inversions in the direction of s, the s-inversion with fixed point xq is uniquely characterized 
by the associated harmonic function h. Thus, the family of inversions we obtain is a one 
parameter family of involutions indexed by the fixed point Xq. In Section 3, we state our 
main result. That gives the path construction of X* in terms of the inverse of X in the 
direction of s with respect to a point Xq G E. Section 4 is devoted to applications. We 
point out in Corollaries [2] and [3] some new results that we obtain for Bessel processes and 
the hyperbolic Bessel process of dimension 3. 

2. Conditioned diffusions and inversions 

Let X G VJ^{E). To start with, assume that X is of type 1 and let h : E ^ R+ be a 
positive harmonic function for the infinitesimal generator L of X which is not constant. 
Let X* be the dual of X with respect to h{x)m{dx). We are ready to state the following 
result which motivates the construction of inversions; to our best knowledge, the described 
role of Xq and Xi for X and X* has not been known. 

Proposition 1. 1) There exists a unique xq & E such that h'^{xo) = h{l)h{r). We call 
Xq the h-geometric mean of {l,r}. Furthermore, for x G E, we have F^iHi < Hr) = 
P*(ifr < Hi) if and only if x = Xq. 

2) There exists a unique Xi E E such that 2h{xi) = h{l) + h{r). We call Xi the h- 
arithmetic mean of {l,r}. Furthermore, for x E E , we have Px-(-f^i < Hr) = Px(-f^r < 
Hi) = 1/2 if and only if x = Xi. 

Proof. 1) Since h is continuous and monotone, because it is an affine function of s, if 
h is increasing (resp. decreasing) then the inequality h{l) < \Jh{l)h{r) < h{r) (resp. 
h{r) < ^h{l)h{r) < h{l)) implies the existence and uniqueness of xq. Because — is a 
scale function for X*, see for example [3], applying f ll.3p yields 

¥^{Hi < Hr) = {h{x) - h{r)) I {h{l) - h{r)) 

and 

K{Hi > Hr) = {l/h{x) - l/h{l)) / {l/h{r) - l/h{l)) . 
These are equal if and only if a; = Xq. 2) The proof is omitted since it is very similar. □ 
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To introduce the mappings we are interested in, let us express r in terms of Xq and /, 
when xo is the /i-geometric mean, to find r = {h'^{xo) /h{l)). This exhibits the function 
/ : X — !■ h^^ {h'^{xo)/h{x)) which is well defined on E by monotonicity of h. Clearly, J is a 
decreasing involution of E with fixed point xq. Next, observe that hoIoh~^ : x — /i^(xo) / x 
is the Euclidian inversion with fixed point h^xo). We return now to the general case and 
assume that X is of either of the types 1-4. Our aim is to determine the set of all 
involutions associated to X which lead to the set of Mobius real involutions 

MI := {lu : M\{a/c} R\{a/c}; uj{x) = {ax + b) /{cx - a), + 6c> 0, a, 6, c G M} . 

Note that the condition + 6c > for u G AiT ensures that u, when restricted to either 
of the intervals (—oo, a/c) and (a/c, +oo), is a decreasing involution. Let us settle the 
following definition. 

Definition 1. Let s and be, respectively, a reference scale function for X and its 
inverse function, and Xq E E. A mapping I : E ^ E is called the inversion in the 
direction of s, or s -inversion, with fixed point Xq if the following hold: 

1) / o /(x) = X for X G E; 

2) solos'^ e MI; 

3) I{E) = E; 

4) I{Xo) = Xq. 

If s o I o is the Euclidian reflection in Xq then I is called the s-reflection in Xq. 

Since I is defined on the whole of the open interval E, it is necessarily continuous. 
This, in turn, implies that it is a decreasing involution such that /(/) = r. The objective 
of our next result is to show the existence of the inversion of E in the direction of s in 
case when s is bounded on E i.e. for diffusions of type 1. 

Proposition 2. Let xq G £" and assume that s is bounded on E. Then, the following 
assertions hold. 

1 ) The inversion of E in the direction of s with fixed point xq is given by 




V 
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where 

a = {2s{l)s{r) - s(xo) (s(0 + s{r))) {s^{xo) - s(/)s(r))-' s(xo) 

b = (2.(xo) - + s(r))) (s2(xo) - s(0^(r))"' • 
2) If2s{xo) ^ + s{r) then I = h~^{l/h) where 

{bs{x) - l)/ibs{xo) - 1) ifs^ixo) ^ s{l)s{ry, 
s{x)/s{xo) otherwise. 

Furthermore, h is continuous, strictly monotonic and satisfies h ^ on E. 



h{x) 



Proof. 1) We shall first assume that s'^{xo) ^ s{l)s{r). We look for / such that solo 
s~^{x) = {x + a)/ [bx — 1) where a and b are reals satisfying a6 + 1 7^ 0. Since the images 
of / and Xq by / are respectively r and Xq, we get the following linear system of equations 



(2.6) 



6s^(xo) — a = 2s(xo); 
bs{l)s{r) — a = s{l) + s{r). 

Solving it yields a and b. We need to show that ab + 1 ^ 0. A manipulation of the first 
equation of (12.61) shows that 1 + ab = {bs{xo) — 1)^. In fact, we even have the stronger fact 
that s{x) 7^ 1/b on E which is seen from 1/6 > s(r) if 2s{xq) > s{l) + s{r) and 1/6 < s{l) 
if 2s{xo) < s{l) + s(r). Finally, if s^(xo) = s(/)s(r) then clearly I{x) = s~^{s'^{xo)/s{x)). 

2) Assume that 2s{xo) 7^ s(/) + s(r). Let us first consider the case s^(xo) 7^ s{l)s{r). 
Setting h{x) = {s{x) — 1/b) / 6 we then obtain 

h-\l/h{x)) = s-^ {{s{x) + {bS^ - l/b))/{bs{x) - 1)) . 

Thus, the equality I{x) = h~^{l/h) holds if and only ii 5 = ±-\/l + ab/b which, in turn, 
implies that h{x) = ±{bs{x) — 1) / {bs{xo) — 1) . Since h is positive, we take the solution with 
plus sign. Since h is an affine transformation of s, it is strictly monotone and continuous 
on E. Finally, because s ^ 1/b, as seen in the proof of 1), we conclude that h does not 
vanish on E. The case s^(xo) = s(/)s(r) is completed by observing that this corresponds 
to letting 6 — 7- 00 and 6 = s{xo) above which gives the desired expression for h. □ 

Now, we are ready to fully describe the set of inversions associated to the four types 
of diffusions described in the introduction. The proof of the following result is omitted, 
keeping in mind that when s is unbounded on E, by approximating E hj a family of 
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intervals (a, (3) d E where s is bounded, using continuity and letting a — I and /3 — )■ r 
we obtain an expression for /. 



Proposition 3. All kind of inversions of E in the direction of s with fixed point xq E E 
are described as follows. 

1) X is of type 1 with 2s{xq) ^ s{l) + s{r) then the inversion is given in Proposition^^ 

2) X is of type 2 then we have 

I{x) = {s{l) + {s{xo) - s{l)f/{s{x) - .(/))) . 

3) X is of type 3 then we have 

I{x) = {s{r) — (s(r) — s(xo))^/(s(r) — . 

4) X is of type 4 or type 1 with 2s{xq) = s{l) + s(r) then I is the s-reflection in Xq. 

Going back to Proposition [1] we can express the inversions of Proposition |3] in terms 
of the harmonic function h instead of the reference scale function s. For that we need to 
compute the positive harmonic function h described in the introduction for each of the 
types 1-4 of diffusions. We easily get 



(2.7) h{x) 



bs{x)- 


-1 


f)s(X())- 


-1 


1 




s{x)~- 


<l) 


s{xo)- 


s{l) 


s{r)—^ 





if X is of type 1 and 2s(xo) 7^ s{l) + s(r); 
if X is of type 1 and 2s(xo) = s{l) + s(r) or type 4; 
if X is of type 2; 
if X is of type 3. 



V s{r)-s{xQ) 

Note that the case where X is of type 1 and xq is the s-geometric mean is covered in the 
first case by letting 6 —t- 00 to obtain h{x) = s{x)/ s{xq). In the following result, which 
generalizes Proposition [H we note that the first assertion could serve as the probabilistic 
definition for s-inversions. 

Proposition 4. 1) A function I : E ^ E is the s-inversion with fixed point xq E E if 
and only if I{E) = E and for all x E E 

(2.8) ¥,,iH, < Hj^,)) = P:„ (///(.) < H,,) 

where P* is the distribution of the Doob h-transform of X by some positive harmonic 
function h. Furthermore, The s-inversion and the h-inversion of E with fixed point Xq 
are equal. 
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2) Let Q^.(, be the probability measure of {I{Xt),t < Q when Xq = Xq. Then, clearly, 
formula 1^2. 8i) holds true when is replaced by Qxq- We call the process {I{Xt),t < Q 
the inverse with respect to xq of {Xt,t < (). 

Proof. 1) If Xq is the s-arithmetic mean of {/, r} or X is of type 4 then we are looking for 
I : E ^ E such that F^^^Hx < Hj^x)) = 1/2. Using (O} we get {s{xo) - so I{x))/{s{x) - 
s o I{x)) = 1/2 which gives that / is the s-reflection. For the other cases, using ( 12.80 
and the fact that —1/h is a scale function for X*, we find that I{x) = {h{xo)'^ / h{x)) 
so that / is an /i-inversion with fixed point Xq. The "only if" part is straightforward 
following a similar reasoning to that of the proof of Proposition [1] giving I to be either 
the /i-refiection or the /i-inversion with fixed point xq. 

2) This is easily seen by using the first assertion. □ 

For completeness, we show in the following Proposition how to define rigorously con- 
ditioning a transient diffusion to hit one boundary of an interval before another. For 
conditioning a transient diffusion to avoid one of the boundaries we refer, for example, to 

m, m)- 



Proposition 5. Assume that X is transient. Let h be as given by {2/1} and X* be the 
dual of X with respect to h{x)m{dx). Let p be the probability that X, when started at Xq, 
hits I before r. X* is X conditioned to hit I before r with probability q = 1 — p. 

Proof. By construction, we have /i(xo) = 1. Assume at first that X is of type 1. Let us 
decompose h, in terms of hi and hr which are defined below, as follows 

7/ N *Hx)-h{l) h(r) - h(x) *, / N *, / N 

''<^' = '' mI„)-m/) +P m.) -mI„) 

where 



and 



_ h{r)-h{xo) _ h*{r)-h*{Xo) ,rr . rr^ 
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But, we have that q* = p and p* = q when X and X* are started at xq. Thus, for any 
bounded J-j^-measurable functional G and t > 0, we can write 

El^[G{X,,s<t),t<C] = E,MXt)G{Xs,s<t),t<C] 

= qE,,[hi{Xt)G{Xs,s<t),t<C] 
+ pE,,[hr{Xt)G{Xs,s<t),t<Q. 

Next, since our assumptions imply that p 

E,,[hi{Xt)G{Xs,s<t),t<C] = 

where we used the strong Markov property for the last equality. Similarly, for the other 
term, since q G (0, 1) we get 

E,,[{hr{Xt)G{Xs,s< t),t < C] = E^^[G{Xs,s< t),t< C\Hr < Hi]. 

The last two equations imply our assertion. Assume now that h{r) = oo. Then h{l) = 
and PxQ-a.s. all trajectories of the process X tend to / and p = Fxg{Hr < Hi) = 0. We 
follow [18] to define X conditioned to avoid / as follows. For any bounded J-'/'^-measurable 
functional G and t > 0, we set 

E:jG(X„s<t),t<C] = \imE,,[G{X,,s<t),t<C\Ha<Hi] 

a—>-r 

= \imE,,[G{Xs,s< t),t< Ha < Hi]/¥,,,{Ha < Hi) 

a—>-r 

= lim E,.„ [hr{Xt)G{X,, s<t),t<HaAHi] 
= E^,[h{Xt)G{Xs,s<t),t<C] 

where we used the strong Markov property for the third equality and the monotone 
convergence theorem for the last one. In this case P*^-a.s all trajectories of the process 
X* tend to r and p* = P*p(/Jz < Hr) = which completes the proof of the statement. 
The case h{l) = — oo and h[r) = can be treated similarly. □ 

Remark 1. Observe that the s-inversions described in Proposition\^ solve G{x,y) = in 
y, where G is the symmetric function G{x, y) = As{x)s{y) — B{s{x) + s{y)) — C for some 
reals A, B and G . This is in agreement with the fact that I is an involution, see [2T] . 



= f^^{Hi < Hr) G (0, 1), we have 
E,JG{X,,s<t),t<C\Hi<Hr] 
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Remark 2. The inversion in the direction of s with fixed point Xq does not depend on the 
particular choice we make of s. Tedious calculations show that the inversion of E in the 
direction of s is invariant under a Mobius transformation of s. 

Remark 3. From the point of view of Martin boundaries, the functions hi and h^. which 
appear in the proof of Proposition are the minimal excessive functions attached to the 
boundary points I and r, see (^[3], [16], That is the Doob h-transformed processes 

using hi and hj. tend, respectively, to I and r. Harmonic functions having a representing 
measure with support not included in the boundary set of E are note considered in this 
paper since we do not allow killings inside E. 



Let X G VJ^{E) and s be a scale function for X. For xq ^ E, let I : E —> E 
be the inversion of E in the direction of s with fixed point Xq. Let h be the positive 
harmonic function specified by (12. 7p . Let X* be the dual of X with respect to h{x)m{dx). 
As aforementioned, the distribution of X* is obtained as a Doob /^-transform of the 
distribution of X by using the harmonic function /i, as given in (11 ■4p . Clearly, if X is of 
type 1 (resp. of type 2 and drifts thus to /, of type 3 and drifts thus to r or of type 4) then 
X* is of type 1 (resp. of type 3 and drifts thus to r, of type 2 and drifts thus to / or of type 
4). Recall that, for a fixed t < (, Tt is the inverse of the strictly increasing and continuous 
additive functional At := £ r^{Xs)a'^{Xs)/a^ o /(X,) ds; r; and A* are the analogue 
objects associated to the dual X*. Recall that the speed measure m{dy) = 2dy/{a'^s') of 
X is uniquely determined by 



for any J = {a, (5) C. E and all x,y & J, where c is a normalization constant and Gj{., .) is 
the potential kernel density relative to m{dy) of X killed when it exits J; see for example 
[H] and [15]. Recall that —1/h is a scale function and h'^{x)m{dx) is the speed measure 
of X*. We are ready to state the main result in this paper. The proof we give is based 
on the resolvent method for the identification of the speed measure, see ([E], [15], [T9]). 
Other possible methods of proof are commented in Remarks H] and [5l 



3. Inversion of diffusions 




where 



Gj(x, y) = c{s{x Ay)- s{a)){s{f3) - s{x V y))/{s{f3) - s{a)) 
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Theorem 1. With the previous setting, let I be the s-inversion of E with fixed point 
xq G E. Assume that Xq, Xq G E are such that I{Xq) = Xq. Then the following 
assertions hold true. 

(1) For all t < C Tt and A^ have the same distribution. 

(2) The processes {X^,t < (*) and {I{X^^),t < A(^) are equivalent. 

(3) The processes {Xt,t < Q and {I{X*,),t < A*^,) are equivalent. 

Proof. 1) Let t > be fixed and set rjt = I{Xrt). Because Tt is the inverse of At, we can 
write At-^ = t. Differentiating and extracting the derivative of Tt yields 

d a^oUXr,) 



dt r\x,,)a\Xr,y 

Integrating yields 



Tt= f {I'a/{aoI))-\XrMs = Al 
Jo 



'0 

The proof of the first assertion is complete once we have shown that rj and X* have the 
same distribution which will be done in the next assertion. 

2) First, assume that h is not constant. In this — l/h{x) is a scale function 

for 7] since — l//io/(X^J = —h{X^^) is a continuous local martingale. Next, let J = [a, /3) 
be an arbitrary subinterval of E. We proceed by identifying the speed measure of t] on J. 
By using the fact that the hitting time H'^ oi y hj t] equals Ahj^^^ for y ^ E, we can write 



dAr 

Gj^j){I{x),y)I'\y)a\y){a^oI{y)Y^m{dy) 
= 2 r'^^ Gj^j){I{x),y)I'\y){a'oI{y)s\y)YUy 

Jl(l3) 

= 2 Gj^j){I{x),Iiy)){a'iy)i-hoiyiy)}'Uy. 

J a 

On the one hand, we readily check that {a'^{y){s o iy(y)}^^dy = h'^{y)m{dy) = m*{dy) 
for y & J. On the other hand, we have 

{h{I{x) A I{y)) - himmiHa)) - h{I{x) V I{y))) 



Gj^j){I{x),I{y)) 



h{I{a)) - h{im 
{-h*{x V y) + h*{(3)){-h*{a) + h*{x A y)) 
-h*{a) + h*{/3) 

G*j{x,y), 



13 

where h* = —1/h and Gj is the potential kernel density of X* relative to rrf^dy). The 
case when h is constant can be dealt with similarly but by working with s instead of h. 
This shows that the speed measure of rj is the same as that of X* in all cases. Now, since 
= we get = which, in turn, implies that Af = . Finally, using the fact 
that I is an involution gives 

H; = M{s > 0, 1{XAt) = y} = r^,^^, = Ah,,^, 

for y E E. The assertion is completed by letting y tend to either of the boundaries to find 
C* = Aq and C = A*^, as desired. 

3) The proof is easy using 1) and 2), the fact that / is an involution and time changes. 

□ 

If s is bounded at one of the boundaries and unbounded at the other, i.e. X is either of 
type 2 or of type 3, then formula (11. 3p gives that one of the boundaries is almost surely 
hit before the other one. The lifetimes of X* and X are in this case given as follows. 

Corollary 1. If X is of type 2 then ( = Hi and Aj^ have the same distribution and 
(* = CO. Similarly, if X is of type 3 then (* = and A^o have the same distribution 
and ( = oo. 

Remark 4. Since X satisfies the s.d.e. ( (i. ij) . by Girsanov's theorem, we see that X* 
satisfies Yt = X* + J^a{Ys)dBs + /o(6 + a^h' /h){Ys)ds for t < C where B is a Brown- 
ian motion which is measurable with respect to the filtration generated by X* and C,* = 
inf{s, Yg = I orYg = r}, see for example [6]. Long calculations show that rj also satisfies 
the above s.d.e. which, by Engelbert- Schmidt condition ( fi.^) . has a unique solution in 
law. This gives a second proof of TheoremUl Note that the use of Ito's formula for rj is 
licit since I G C'^{E). 

Remark 5. Another way to view the main statements of Theorem U\ is the equality of 
generators 

^ ""'^"^^ L(/o/)(/(x)) = -l-W)(x) = L7(x) 



J'2(x)a2(J(x)) ' h{x 

and 

for all X E E, g E T^i^X) and f G However, the main difficulty of this method of 

proof is the precise description of domains of generators. 
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4. Applications 

4.1. Inversions of Brownian motions killed upon exiting intervals. Assume that 
A is a Brownian motion killed upon exiting the interval E = {l,r). Let Aq = Xq ^ E. If 
E is bounded then we obtain the inversions 

xl/x, if xl = Ir; 

2xq — X, if 2xo = / + r; 

(x + a)/{bx — 1), otherwise, 



Kx) 



where 



2/r-Xo(/ + r) ^ , 2xo - {I + r) 

-Xo and o — 



Xq — Ir Xq — Ir 

where i? is a standard Brownian motion. We distinguish three cases appearing in the form 

of the inversion /. In the first case equation x^ = Ir implies that I and r are of the same 
sign. If / > (resp. / < 0) then A* is the (resp. negative) three dimensional Bessel process 
killed upon exiting E. In the second case, A* is a Brownian motion killed when it exits E. 
In the third case. A* satisfies the s.d.e. = Bt+xo + J^{X* — l/b)~^ds for t < (*. If xq is 
below the arithmetic mean i.e. Xq < (/+r)/2 then xq — 1/& = {xo — r){xo — l)/{2xo — {l+r)). 
By uniqueness of the solution to the s.d.e. Rt = Ro + Bt + J^{l/Rs)ds deriving the three- 
dimensional Bessel process R, we get that A^* = l/h + Rt with Rq = xq — 1/h killed when 
R exits the interval (/ — 1/6, r — 1/6). If xq is above the arithmetic mean then we find 
XI = 1/6 — Rti where i? is a three-dimensional Bessel process started at 1/6 — xq, where 
R is killed as soon as it exits the interval (1/6 — r, 1/6 — /). 

If / is finite and r = oo or / = — oo and r is finite then by Proposition [3] we respectively 
obtain 

/(^) = / + 1^11^ and I{x) = r - t^i^. 

X — I r — X 

If r = +00 then a similar reasoning as above gives that Aj* = I + Rt where i? is a three- 
dimensional Bessel process started at xq — I. If / = — oo we obtain X^ = r — Rt where 
i? is a three-dimensional Bessel process started at r — xq. If 2xo = / -|- r or / = — oo and 
r = -|-oo then we obtain the Euclidian refiection in xq, i.e. x — 2xo — x, and A* is a 
Brownian motion killed when it exits E. Note that for E = (0, oo) the conclusion from 
our Theorem[T]is found in Lemma 3.12, p. 257, of |13]. That is (A^, t < Hq) is distributed 
as (1/A**, t < A*^), where r* is the inverse of = J*{X*)~'^ds. In this case A* is the 
three dimensional Bessel process; see our last Example 4.3 for Bessel processes of other 
dimensions. If ii^ = R then / is the Euclidian refiection in xq. Finally, observe that the set 
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of inversions of E we obtain for the Brownian motion killed when it exits E is precisely 
MI{E) := {/ G MI; I{E) = E}. 



4.2. Inversions of drifted Brownian motion and hyperbolic Bessel process of 
dimension 3. Set bI-^^ = Bt + jit, t > 0, where i? is a standard Brownian motion and 



if /i > (resp. /i < 0). Let us take the reference scale function s{x) = — e~^^^/(2yu). 
Observe that s is increasing for all fi ^ 0. Moreover lima,_j.oo ^(x) = if /i > and 
lima;_>._oo s{x) = if /i < 0. We take X to be B^^^ killed when it exits (/,r) C M. Let us 
fix xo G E. 

If we take E = M. then by Proposition |3l even though X is of type 2 if /i < and of 
type 3 if /i > 0, the inversion of E in the direction of s is the Euclidian reflection in xq. 
X* is the Brownian motion with drift fi* = —fi in this case. If s{l) and s(r) are finite then 
using Proposition |2] we obtain appropriate, but in most cases complicated, formulas for 
I. For Brownian motion with drift the case of the half-line is the most interesting. Take 
for instance E = (0, oo) and process Xt = Xq + Bt + fit starting from some point Xq > 
and killed at zero. We consider two cases: if /i < then s(0) = —l/{2fi) > 0, s(oo) = oo 
and X is of type 2; if > then s(0) = —1/(2//) < 0, s(oo) = and X is of type 1. 

First, let fi < 0. If X is of type 2 then, by Proposition |3l we have only one possible 
inversion: /(x) = + (s(xo) — s(/))^/(s(x) — which gives 



Now, if /i > then X is of type 1. Because e"^'^^ 7^ implies s^(xo) 7^ s(0)s(oo) = only 
two cases are possible. Either 2s(xo) = s(0) + s(oo) = s(0), which gives Xq = 2^ In 2, and 
then we have s-reflection J(x) = —j^ ln(l — e"^'^^) or 2s(xo) 7^ s(0) + s(oo) and then the 
formula from Proposition 3 gives the inversion 



/i G M, /i ^ 0. Thus, 



is a transient diffusion which drifts to +00 (resp. to —00) 



/(x) = (2|//|)-Mn ((e-2^" - 1 + (1 - e^2'^"°)')/(e""'" - 1)) • 



If we choose Xq = (2|/i|) ^ln(l + -\/2) then the above formula simplifies to 



J(x) = (2|/i|)-Mn ((e"''^" + l)/(e-2^" - 1)) = (2|/i|)-Mn coth(|/i|x). 




This in turn simplifies if we choose e^^^° — 2e^'^^'° = 1, that is, if xq 



ln(l + y2). Then 
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It is easy to check that if yU < 0, then h{x) = ^ a-nd then X*, being an /i-process, 

has generator L* f{x) = + coth(|/_i|x)/'(x). In particular, if /i = — 1 then 

I{x) = ^Incothx and X* has generator L*f{x) = + coth(x)/'(x). This recovers 

the well-known fact that Bt — t conditioned to avoid zero is a three-dimensional hyperbolic 
Bessel process. The novelty here is that we get X* as a time changed inversion of Bt — t. If 
/i > 0, then h{x) = ^ ^^^"^ and X* has the generator L*f{x) = + yutanh(/ix)/'(x). 

Note that we have shown the following particular result. 

Corollary 2. Let X be a three-dimensional hyperbolic Bessel process started at Xq = 
iln(l + v^) and I{x) = |lncothx. Then, I{X) is a time-changed drifted Brownian 
motion Bt — t conditioned to avoid zero. In particular, the functional 



POO 

Aoo = (1/4) / ds/ (cosh(X,) sinh(X,))^ ds 
Jo 



has the same distribution as the first hitting time of by the Brownian motion with minus 
unity drift. 

4.3. Inversions of Bessel processes. It is known that the Bessel processes of dimen- 
sions 6 and 4 — 6 are dual one to another, see e.g. ([3], [6], [H]). But, except the case 
5 = 1, it was not known that, up to a time change, they can be obtained one from another 
by the deterministic inversion x — > 1/x. 

Let X be a Bessel process of dimension 6 >2. Thus, E = (0, oo), is polar and X has 
the infinitesimal generator Lf{x) = \f"{x) + {5 — l)/(2x)/'(x) for x E E, see ([H], [E])- 
Let z/ = (5/2) — 1 be the index of X. The scale function of X may then be chosen to be 

-x-"^" if 1/ > 0; 
s[x] = < 21oga; if = 0; 

x~'^'^ otherwise. 

X is recurrent if and only if 5 = 2, see for example Proposition 5.22, p. 355 in [9]. First, 
when 6 = 2 the inversion of E in the direction of s is the s-refiection x — )■ Xq/x. For 
the other cases, the inversion of E in the direction of s with fixed point 1 is found to be 
a; — 1/x. Furthermore, X* is a Bessel process of dimension 4 — 6 < 2. Observe that in 
the two considered cases the involved clock is At = /^(xo/Xs)^ ds. Hence, we have shown 
the following result. 
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Corollary 3. Let X be a Bessel process of dimension 5 G M, killed at if 6 < 2, starting 
from Xq > and I{x) = 1/x. Then I{X) is a time-changed Bessel process of dimension 
4 — 6. In particular, the functional ds/{XsY, when the dimension of X is larger than 
1, has the same distribution as the first hitting of zero by the Bessel process of dimension 
4-6. 

Acknowledgment: We are greatly indebted to I'Agence Nationale de la Recherche for 
the research grant ANR-09-Blan-0084-01. 
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